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ABSTRACT 

Navier-Stokes simulations of high Reynolds num- 
ber flow around an axisymmetric body supported in 
a water tunnel were made. The numerical method is 
based on a finite- differencing high resolution second- 
order accurate implicit upwind scheme. Four different 
configurations were investigated, these are, bare- 
body, 21 body with an operating propeller, 3.) body 
with a ring wing and 4. body with a ring wing and 
an operating propeller. Pressure and velocity compo- 
nents near the stern region were obtained computa- 
tionally and are shown to compare favorably with the 
experimental data. The method correctly predicts the 
existence and extent of stern flow separation for the 
barebody and the absence of flow separation for the 
three other configurations with ring wing and/or pro- 
peller. 

I. INTRODUCTION 

It has been known in marine propulsion technol- 
ogy that certain advantages can be achieved by prop- 
erly integrating a ring wing and the propeller. With 
a ring wing that accelerates the flow the efficiency of 
the propeller remains more or less at a constant level 
over a wider range of advance ratios. With a ring 
wing that decelerates the flow, the inception of cav- 
itation on the propeller can be delayed. In order to 
take advantage of the ring wing to obtain the desir- 
able benefit, a clear understanding of the role it plays 
is important. Thorough water tunnel experiments and 
reliable numerical simulations serve as complementary 
approaches to gain understanding. 

By combining the body, ring wing and propeller, 
appropriate configurations can be generated for water 
tunnel experiments, and the influence of each individ- 
ual component can be isolated. The configurations 
studied here are 1. barebody, 2. body with an oper- 
ating propeller, 3. body with a ring wing and 4. body 
with a ring wing and an operating propeller. Data 
collected from the water tunnel tests included the ve- 
locity components around the afterbody and the pres- 
sure on the stern. The data show the degree of inter- 


action at the given operating conditions and serve as 
benchmarks for evaluating the present numerical sim- 
ulations. 

One important parameter in hull-propulsor inter- 
action is the thrust deduction coefficient t, which signi- 
fies the drag augmentation due to the interaction. In 
the past, inviscid methods [1,2] have been successful 
in computing the coefficient t. The methods become 
somewhat inadequate in a situation where the propul- 
sor unit is imbedded in the stern boundary layer where 
the viscous effect plays a dominating role. Efforts 
have been made to address the problem; in particular, 
Falcao de Campos [3] presented an inviscid approach 
to calculate the flow on the stern with and without 
propulsor based on the Euler equation of motion and 
Huang et al [4] developed a numerical technique to 
study the interaction between a propeller and unsep- 
arated viscous stern boundary layer. To further en- 
hance the ability to predict the effect of hull-propulsor 
interaction, a Navier-Stokes type viscous analysis is 
needed. This is particular true when barebody flow 
separation may occur. Previously, Haussling et al [5] 
performed extensive numerical simulations of viscous 
flow about bodies with appendages using a Navier- 
Stokes solver. Here a three-dimensional incompress- 
ible Navier-Stokes solver is used to simulate the flow 
around a compound propulsor unit on an axisymmet- 
ric body supported in a water tunnel with a square 
cross-section. The solver is based on a high resolution 
second-order accurate implicit upwind scheme [6,7]. 
The propeller effect is simulated by imbedding body 
forces in a disk located at the propeller plane [8,9]. The 
experimental data were used to validate the Navier- 
Stokes solver. 

II. DESCRIPTION OF EXPERIMENTS 

The test body is axisymmetric with a length of 
139.12 cm and a maximum diameter of 24 cm; Its 
radius offsets, y, as a function of axial length, x, are 
given in cm units by: 
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for 0.00 < x < 24.00 

y = [12.0 2 - (12.0 — x/2) 2 ] 1 / 2 
for 24.00 < a: < 97.81 
3/ = 12 

for 97.81 < x < 115.84 

y = [42. 67 2 - (s - 97.81) 2 ] 1 / 2 - 30.67 
for 115.84 < x < 129.17 
y = (133.0 - x) tan25° 
for 129.17 < x < 131.93 
y = —0.0214a; -f 4.554 
for 131.93 < x < 139.12 

y - -0.0310a; 2 + 8.162a; - 535.51 

The ring wing has a NACA 4415 profile section 
with a 5° angle of attack; its chord length is 5.3 cm. 
The diameter of ring wing measured at its trailing 
edge is 16 cm; its trailing edge is located at 129.7 cm 
from the nose of the body. The propeller has four 
blades with a diameter of 15.71 cm and was driven 
from behind with a Z-drive propelling device. The de- 
vice consists of a tapered forebody, a cylindrical mid- 
body and an elliptical afterbody. The length of the 
device is about 85 cm and the diameter of midbody 
is about 10 cm. The propeller plane is located at 135 
cm from the nose. The water tunnel in which the ex- 
periments were conducted has a square cross-section 
with round corners and its dimensions are 90 cm X 90 
cm. The model/tunnel blockage ratio is about 5.6 per- 
cent. The body was supported from the ceiling of the 
water tunnel with two struts located at x = 14.0 cm 
and x — 64.0 cm from the nose. A schematic sketch is 
shown in Fig. 1.; the propeller drive is not included. 
The geometric tolerances are less than 0.5 mm. The 
reference length is chosen to be 133.0 cm. 

Experiments were conducted at a Reynolds num- 
ber of 6 x 10 6 (based on reference length 133.0 cm). 
The flow measurements were carried out on lines lying 
on the horizontal plane and perpendicular to the cen- 
ter line. The pressures on the body surface were mea- 
sured by means of transducers and all signals could 
be reproduced satisfactorily within 1.0 percent. For 
the flow measurements, the standard deviation of all 
signals varied between 0.1 and 2.0 percent. 

III. DESCRIPTION OF NUMERICAL 
PROCEDURE 

Using Chorin’s artificial compressibility formula- 
tion, the incompressible Navier- Stokes equations are 
written in conservation law form for three-dimensional 
flow as [10] 

Q< + (E* - e;) x + (F* - f;) 9 + (G* - G*)* = 0 (1) 

where the dependent variable vector 

Q = (p,«,v,w) T 


represents the pressure and velocity components in a 
Cartesian coordinate system (a;, y,z), The inviscid flux 
vectorsE*, F*, G* and the viscous shear flux vectors 
E* , F* , G* are given by 

E* = (fiu,u 2 + p,uv,uw) T 
F* = {fiv,uv,v 2 +p, vw) T 
G* = (0w,uw,vw,w 2 -f p) T 
E v — Re (0, Txx5 T xyi T xz) 

= Re 1 ( 0 , Ty x , Tyy , Ty z ) ^ 

G* = Re” 1 (0, Tzxi r zy , t zz ) t 

The coordinates x,y>z are scaled with an appropriate 
characteristic length scale L. The velocity components 
u,v y w are nondimensionalized with respect to the free 
stream velocity while the normalized pressure is 
defined as p = (P — Poo)/ pVoo- The kinematic vis- 
cosity v is assumed to be constant, and the Reynolds 
number is defined as Re — The artificial com- 

pressibility parameter 0 monitors the error associated 
with the addition of the unsteady pressure term in 
the continuity equation which is needed for coupling 
the mass and momentum equations in order to make 
the system hyperbolic. 

Equations (1) can be transferred to a curvilinear, 
body-fitted coordinates system (£ , £ and rj ) through 
a coordinate transformation of the form 


C = COv?/, z), i = ^x,y,z) and r) = »?(*, y,z) 
Eq. (1) becomes 


(Q/J)t + (E - E.) { + (F - F v )< + (G - G,), = 0 (3) 


with 

and 


, where 


(E,F,G) t = [T] (E*,F*,G*) T 
(E, F, G) t = [T] (E*,F*,G*) r 


[T] = 


G 

G 

rjx 


Cy G 

G G 

Vy 11 z 


and the Jacobian of the coordinate transformation is 
given by 



y< 


x € 

vs 

H 

Xy 

y n 

z v 


The Cartesian derivatives of the shear fluxes are ob- 
tained by expanding them using chain rule expansions 
in the £, £, and r\ directions. 


Defining computational cells with their centroids 
a t / — ^ ( 9 is £, £,or rj) and their cell interfaces 
at l ± 1/2, the backward Euler time differencing of the 
three-dimensional conservation form is 
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- == - [A c (E n+1 - E" +1 ) + A*(F n+1 - F" +1 ) 

+ A,(G " +1 - G” +1 )] 

( 3 ) 

where At is the time step, AQ n = Q n+1 — Q u and 
A/( ) = [( )i+i/2 — ( )i- 1/2]/ AO . Superscript denotes 
the time level at which the variables are evaluated. 

Linearizing Eq. (3) about time level n, we obtain 


ation is used in the streamwise direction ( in combina- 
tion with approximate factorization in the remaining 
two coordinate directions f and rj. It is used to avoid 
the A t 3 spatial splitting error incurred in fully three- 
dimensional approximate factorization methods. This 
scheme is unconditionally stable for linear systems and 
offers the advantage of being completely vectorizable 
like a conventional three-dimensional approximate fac- 
torization algorithm. As a result, Eq. ( 5 ) becomes 


[JL + { m 

l A tJ K dQ 


dE n dF n 


9G» 

w )A( 


+( 


dG n 


dG n 

f )A,]AQ B 


dQ dQ 

= _ [A C (E» - FJ) ■+ . A*(F B - FJ) + A„(G" - G£)] 

( 4 ) 

where I is the identity matrix. 


The left hand side is the implicit part and the 
right hand side is the explicit part of the formulation. 
The explicit part is the spatial derivatives in Eq. 2 
evaluated at the known time level n ; its value di- 
minishes as the steady state solution is approached. 
Hence , it is also called the residual. The L 2 norm of 
the residual is often used as a measure of convergence 
of a solution. Letting the flux Jacobians A, B and C 
be defined as follow 


A = 


dE n __ dF n _ dG n 

<9Q ’ B = dQ C ~ dQ ’ 


discretize the inviscid and viscous fluxes according to 
upwind differencing scheme and central differencing 
scheme respectively in (, £ and 77 coordinate direc- 
tion independently and then assemble them together. 
Equation (4) becomes 


-(A" + X) j+ iA i+ i + (A+ + X) i _ i A i _i 
-(B- +Y) j+i A j+ i + (B+tY) H A i4 
-(C- + Z) k+ xA k+k + (C+ + Z) i _iA fc _ i ]"AQ” 

= -RES(Q n ) 

( 5 ) 

where i,j, and k are spatial indices associated with 
the f, rj and C coordinate direction. A ± ,B :t and 
C± are flux matrices split from the flux Jacobians 
A,B and C according to the signs of their eigenval- 
ues. The residual RES{ Q n ) is evaluated with a TVD 
technique together with Roe’s [11] flux- difference split- 
ting scheme, the discretization is third-order accurate. 
Conventional second-order central differencing is ap- 
plied to obtain the viscous flux matrices X, Y and 
Z. Equation(5) is solved by an implicit hybrid algo- 
rithm, where a symmetric planar Gauss-Seidel relax- 


+ Y) J+ i A J+ i + (B+ + A M ]AQ 
= -RES{ Q n ,Q" +1 ) 

[M-(C- + Z) k+ xA k+k + (C+ + Z) fc _iA fc _i]AQ n 

=MAQ 

+1 ( 6 ) 

qn+i = qn + A q » 

with 

M = [ ^7 + (A “ + X)i+i + (A + + X),..*]. 

and the residual on the RHS indicates the nonlinear 
updating of the residual by using Q n+1 whenever it be- 
comes availabe while sweeping in the ( direction back 
and forth through the computational domain. 

For laminar flow computations the coefficient of 
molecular viscosity ji = fn is obtained from Suther- 
land’s law. Turbulence is simulated using the Baldwin- 
Lomax algebraic turbulence model[12]. For turbulent 
flow laminar viscosity coefficients are replaced by 

i* = m +/** 

The turbulent viscosity coefficient fi t is computed by 
using the isotropic, two-layer Cebeci type algebraic 
eddy-viscosity model as reported by Baldwin- Lomax. 
Modifications proposed by Degani and Schiff [13] and 
Hartwich and Hull [14] were implemented. 

IV. DESCRIPTION of body force 
PROPELLER MODEL 

The principle of the body force model is to in- 
troduce the body force terms into the Navier- Stokes 
equation to include the effects of the propeller. The es- 
sential parameters that define the propeller effects are 
the thrust coefficient CV, the torque coefficient Cq, the 
advance coefficient J and the radial circulation distri- 
bution G(r). The same parameters were used to define 
the body force for the propeller model. The thrust and 
torque coefficients are defined as follow: 

T Q 

Ct = Cq = 

where T and Q are thrust and torque, respectively. D 
is the diameter of the propeller. The axial and circum- 
ferential body force per unit volume are obtained from 
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the following equations: 

rt C T R 2 p G(r) 

jb x — p 

4AX fg G(r)rdr 

A. — Ca T (r) 

2rAX j* p h G(r)rdr 

where fb x and fb$ are the body forces per unit 
volume in the axial and circumferential directions, re- 
spectively, Rk and R p are the radii of propeller hub 
and blade tip, respectively, and AX is the thickness 
of the disk. The computed body forces are then incor- 
porated into the right hand side of Eq. 5 and form a 
part of the residual. Only a slight modification to the 
flow solver is needed to accommodate the body force 
type propeller model and there is no need for special 
gridding. 

In reality, the blade circulation distribution G(r) 
depends upon the inflow at the propeller plane which 
in turn is influenced by the blade circulation. This mu- 
tual dependency implies that the body forces fb x and 
fb q which are functions of G(r) should be obtained 
by an iterative procedure. To complete this procedure, 
knowledge of propeller- induced axial and tangential 
velocities, u a and u t is needed. A propeller program 
based on the vortex-lattice lifting- surface method de- 
veloped by Greeley and Kerwin [15] can be used for 
this purpose. The iterative procedure can be described 
as follow: 

1. Calculate the nominal inflow with the Navier- 
Stokes solver. 

2. Obtain the circulation distribution G(r), the 
induced- velocities u a and the thrust and torque 
coefficients Ct and Cq by using the calculated nominal 
inflow as input to the propeller program [15]. 

3. Compute the body forces fb x and fbg 
by using the calculated circulation distribution G(r), 
thrust and torque coefficients Ct and 

4. Obtain the total velocities at the propeller 
plane by using the Navier-Stokes solver with the infor- 
mation obtained in step 3. 

5. Compute the effective wake by subtracting 
the propeller-induced velocities obtained in step 2 from 
the total velocities obtained in step 4. 

6. Obtain an updated circulation distribution, 
propeller-induced velocities , and thrust and torque 
coefficients by using the newly computed effective wake 
as input to the propeller program. 

7. Repeat the process from step 1 to step 6 until 
the total velocities, the body forces and the propeller- 
induced velocities are unchanged. 

It has been shown that this procedure converged 
after two iterations [8], For the purpose of illustration, 
the results presented in following are obtained by using 
measured thrust and torque without any iteration. In 
spite of its simplicity, it was able to predict the flow 


pattern around the propeller disk and the stern region 
rather accurately [8,9]. 

V. GRID generation: 

Based on the configuration shown in Fig. 1, which 
models the experimental setup described in Section II 
above, a 180° sector of the tunnel needs to be mod- 
elled in order to resolve the effect on flow due to the 
supporting struts. This requires a large amount of 
grid points and extensive computational resources. Af- 
ter one computation for the barebody configuration, 
it was found that the struts produced an influence 
around and directly behind them with no significant 
effect on the horizontal plane on which measurements 
were made. Therefore, the struts were eliminated ah 
lowing the numerical simulations to be performed ac- 
curately in a 90° sector of the tunnel. Also included, 
due to its proximity to the stern, is the propeller drive. 

To include the ring wing geometry, two block C — 
O type grids were generated. The grid points were 
matched at the branch cut that separated the two 
blocks. Computations were performed on a CRAY- 
YMP machine which has eight processors. Compu- 
tational efforts on the two blocks can be performed 
on two processors simultaneously. For the purpose of 
synchronization between the processors, it is more ef- 
ficient if the number of computations is the same for 
each block. For this reason, each block has the same 
number of grid points. The grid system was generated 
by a transfinite interpolation technique. Several grid 
systems with different number of grid and distributions 
had been generated and were used for computations on 
barebody configuration in order to investigate the re- 
lationship between the convergence and grid density 
and distribution. The grid systems examined include 
1. a two block 25 X 49 x 110 (r, <9,z) grids, 2. a 
two block 25 X 25 X 90 grids and 3. a two block 
25 x 13 X 98 grids. Grid points were clustered near 
the boundaries such as tunnel wall and body and ring 
wing surfaces where the viscous effect dominates. The 
minimum spacing normal to the body surface for the 
three grid systems mentioned above is 5.0 xlO” 4 . The 
differences between the solutions based on the grid sys- 
tem 1 and 3 are about two percent which is within the 
reported experimental accuracy. The results presented 
below are based on a two block 25 X 13 x 98 (r,0,x) 
coarse grid system. Convergence is achieved when the 
L >2 normal of all residuals is reduced by three order of 
magnitude with CFL=10. Computational CPU time 
is 40 //sec per node per iteration. For each compu- 
tation, over 500 iterations were carried out to ensure 
convergence. 

VI. BOUNDARY CONDITIONS 

On the solid boundaries such as surfaces of tunnel 
wall, ring wing, body and propeller drive the no-slip 
condition is applied, in addition, the normal gradient 
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of the pressure is assumed to be zero. Freest ream con- 
dition is applied as inflow condition. Zero-order ex- 
trapolation is used to obtain the outflow conditions. 

During the computation process, the variables on 
the branch cut that separates the two blocks are not 
computated. At the end of each iteration, these vari- 
ables are updated by averaging the values at the adja- 
cent grid points from each block. The averaging pro- 
cess is linear. The updated values are then used as 
boundary conditions for both blocks for the next iter- 
ation. 

VII. RESULTS AMD DISCUSSIONS 

The four configurations that were investigated are 
1. flow over the barebody, 2. flow over the body with 
an operating propeller, 3. flow over the body with a 
ring wing and 4. flow over the body with a ring wing 
and an operating propeller. Computational results are 
presented in the form of velocity profiles and pressure 
contours in the stern region; they are compared with 
available experimental data. In addition, to facilitate 
the flow visualization and discussion, computed par- 
ticle trace are also included. The length scale was 
normalized with reference length L = 133cm. 

Case 1. Flow over the Barebody 

In this test case, the configuration is simple but 
the flow is rather interesting. At a Reynolds num- 
ber of 6.0 x 10 6 , flow separation was observed in the 
stern region. Figure 2 shows the velocity vectors in the 
stern region from both computation and experiment. 
The correlation is good except at the axial location 
XjL = 0.89 where the experiment shows a some- 
what fuller profile near the body. Figure 3 shows the 
particle traces and clearly depicts the separation bub- 
ble. The predicted separation location is about 1.5 cm 
(1.13 percent of body length) ahead of where it was 
observed experimentally. Figure 3 also indicates that 
the size of the bubble is predicted correctly. Figure 4 
shows the predicted and measured pressure distribu- 
tions on the surface of the stern. The surface pressure 
begins to recover as the flow passes the shoulder of the 
afterbody at XjL — 0.78 . The recovery levels off 
where flow separation takes place. Figure 5 shows the 
pressure contours in the stern region. 

Case 2. Flow over the Body with 
an Operating Propeller 

In this test case, the propeller was operating at 
J = 0.47 and V = 4 m/s; Ct and Cq were measured as 
2.052 and 0.247, respectively. To apply the body force 
propeller model, the circulation G(r) was assumed to 
be distributed over the disk according to: 

G(r) = r(l - r) 1/2 


with r == [y — yhub)/(R P — Vhub)- The propeller was 
located at X/L = 1.015. 

Figure 6 shows the comparison of predicted and 
measured velocity vectors in the stem region. The pre- 
diction confirms that the separation bubble is removed 
due to the propeller suction. Experimental data in- 
dicate stronger downward radial velocities at all four 
stations. In addition, at axial locations XjL — 0.89 
and Xj L = 0.93 the predicted velocity profiles are less 
full than those that were measured near the body. It is 
found from computations that the predictions behind 
the propeller are sensitive to the circulation distribu- 
tion over the propeller disk and the axial locations at 
which the predictions are made. Figure 7 shows the 
particle traces, indicating that flow contraction takes 
place immediately in front of and behind the propeller 
disk. Figure 8 shows the predicted and measured pres- 
sure distributions on the surface of the stern. The dif- 
ference between the pressures presented in Fig. 8 and 
Fig. 4 gives the amount of the pressure drag on the 
barebody due to propeller action. The added drag con- 
stitutes the major part of the thrust deduction fraction 
(t). The agreement between prediction and measure- 
ment is very good. This is an indication that thrust 
deduction (1 — t) can be predicted correctly with this 
numerical method. Figure 9 shows the pressure con- 
tours in the stern region; the pressure jump across the 
propeller plane is clear. 

Case 3. Flow over the Body with a Ring Wing 

An intended function of the ring wing was to ac- 
celerate the flow and produce thrust. It was also ex- 
pected that the flow separation over the stem would 
be removed once the ring wing was in place. Figure 10 
shows the velocity vectors in the stern region from both 
computation and experiment. In comparison with Fig. 
2, it can be seen that the ring wing achieved its func- 
tion in accelerating the flow and removing the stern 
separation. The predicted and the measured velocity 
vectors are in good agreement. Due to lack of details, 
the experimental data failed to resolve the wake struc- 
ture behind the ring wing. Figure 11 shows the parti- 
cle traces and there is no detectable flow separation on 
the surface of either the body or the ring wing. Figure 
12 shows the pressure on the stern surface. Note that 
the predicted surface pressure immediately upstream 
of the ring wing’s leading edge at X/L — 0.93 is some- 
what higher than that measured. Figure 13 shows the 
pressure contours in the stern region. Clearly shown 
are the locations of the pressure and suction peaks on 
the ring wing surface. 
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Case 4. Flow over the Body with a Ring Wing 
and an Operating Propeller 

In this test case, the propeller was operating at 
J = 0.47 and V = 4 m/s; Ct and Cq were measured 
as 2.081 and 0.250, respectively, and are about 2 per- 
cent higher than those shown in Case 2. Figure 14 
shows the comparison of predicted and measured ve- 
locity vectors in the stern region. The agreement is 
quite good. Figure 15 shows the particle traces. It 
can be seen that flow is accelerated as it passes the 
ring wing and is contracted as it passes the propeller 
disk. Figure 16 shows the pressure on the stern sur- 
face. The discrepancy between the predicted and mea- 
sured pressure immediately upstream of the ring wing 
at XjL = 0.93 as discussed in Case 3 is present here 
also. The difference between the pressures presented 
in Fig. 16 and Fig. 12 gives the amount of the pres- 
sure drag on the body due to the propeller action, the 
agreement between the prediction and measurement is 
very good, despite the discrepancy mentioned above. 
Figure 17 shows the pressure contours in the stern re- 
gion. 

VIII. CONCLUSION 

A numerical method based on a finite-differencing 
high resolution second-order accurate implicit upwind 
scheme was used to discretize the three-dimensional 
Navier-Stokes equations. Simulations of flow over an 
axisymmetric body with a compound propulsor sup- 
ported in a square water tunnel were performed. A 
body force type propeller model was used to simulate 
the propeller action. Results indicate that this nu- 
merical method is effective and accurate. Observed 
flow phenomena such as separation, acceleration and 
contraction were realistically predicted. The pressure 
drags due to propeller action were computed correctly. 
A ring wing may affect the circulation distribution on 
the propeller, depending on its proximity to the pro- 
peller. In order to improve the predictions further, 
an iterative procedure described in section IV will be 
explored in a future study. 
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Figure 1. Schematic Sketch of Test Configuration 



Figure 2. Velocity Vectors in Stem Region for Case 1 
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Figure 3. Particle Traces for Case 1 



Figure 4. Pressure Distribution on Stem Surface 
for Case 1 



Figure 5. Pressure contours in Stem Region 
for Case 1 
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Figure 11. Particle Traces for Case 3 


Figure 8. Pressure Distribution on Stem Surface 
for Case 2 
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Figure 12. Pressure Distribution on Stern Surface 
for Case 3 


Figure 15. Particle Traces for Case 4 
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Figure 13. Pressure contours in Stern Region 
for Case 3 


Figure 16. Pressure Distribution on Stern Surface 
for Case 4 
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Figure 14. Velocity Vectors at Stern Region 
for Case 4 
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Figure 17. Pressure contours in Stem Region 
for Case 4 







DISCUSSION 


Philippe Genoux 

Bassin d’Essais des Car&nes, France 

It is well known that Navier Stokes solvers may have difficulties at 
high Reynolds numbers. How does your code behave at full scale 
Reynolds numbers? 

AUTHORS’ REPLY 

The differencing scheme presented in this paper is based on 
hyperbolic formulation, assuming that at high Reynolds numbers the 
behavior of Navier-Stokes equations becomes hyperbolic-like. At 
present, the solver has been used to simulate flow at Reynolds 
number 15 X 10 6 . The agreement between computation and 
experiment is good. No convergence problem was experienced. 


DISCUSSION 
Fred Stem 

The University of Iowa, USA 

The results presented display some interesting features for a ducted 
propulsor; however, it would be of greater interest if the authors 
would include the propulsor-hull interaction, i.e., make full use of 
the viscous-flow approach to propulsor-hull interaction and 
demonstrated in Reference [8], 

AUTHORS’ REPLY 

We appreciate Prof. Stem’s comments and suggestions. We are using 
the method described in Reference [8] for our future work. 
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